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Abstract 

Molecular collisions can be studied at very low relative kinetic energies, in the milliKelvin range, by 
merging codirectional beams with much higher translational energies, extending even to the kiloKelvin 
range, provided that the beam speeds can be closely matched. This technique provides far more intensity 
and wider chemical scope than methods that require slowing both collision partners. Previously, at far 
higher energies, merged beams have been widely used with ions and/or neutrals formed by charge transfer. 
Here we assess for neutral, thermal molecular beams the range and resolution of collision energy that now 
appears attainable, determined chiefly by velocity spreads within the merged beams. Our treatment deals 
both with velocity distributions familiar for molecular beams formed by effusion or supersonic expansion, 
and an unorthodox variant produced by a rotating supersonic source capable of scanning the lab beam 
velocity over a wide range. 
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I. INTRODUCTION 



The frontier field of cold (< 1 K) and ultracold (< 1 mK) gas-phase molecular physics has 
brought forth many innovations [JjJ— . Among motivating challenges is the prospect of study- 
ing collision processes, especially chemical reactions, under "matterwave" conditions. Reaching 
that realm, where quantum phenomena become much more prominent than in ordinary "warm" 
collisions, requires attaining relative velocities so low that the deBroglie wavelength becomes 
comparable to or longer than the size of the collision partners. That has been achieved recently 
for reactions of alkali atoms with alkali dimers formed from ultracold trapped alkali atoms by 
photoassociation or Feshbach resonances fl5l |6fl. With the aim of widening the chemical scope, 
much effort has been devoted to developing means to slow and cool preexisting molecules. (Com- 
pilations are given in flvl— 1 1 Oil . For chemical reactions, however, as yet it has not proved feasible 
to obtain sufficient yields at very low collision energies, using either trapped reactants or crossed 
molecular beams. The major handicap in such experiments is that both reactants must contribute 
adequate flux with very low translational energy. 

Merged codirectional beams with closely matched velocities offer a way to obtain far higher 
intensity at very low relative collision energies, since then neither reactant needs to be particularly 
slow. Moreover, many molecular species not amenable for slowing techniques become available 
as reactants. Merged beams have been extensively used with ions and/or neutrals formed by charge 
transfer, to perform experiments at relative energies below 1 eV with beams having keV energies 
ill ill . A key advantage is a kinematic feature that deamplifies contributions to the relative energy by 
velocity spreads in the parent beams. By virtue of precise control feasible with ions, the velocity 
spreads are also quite small, typically ~0.1% or less. For thermal molecular beams, such as we 
consider here, the spreads are usually ~10% or more. That enables attaining low relative collision 
energy, but much lower energies and improved resolution can be obtained by narrowing the spreads 
to ~1%, which now appears feasible at an acceptable cost in intensity. Surprisingly, application 
of merged beams to low-energy collisions of neutral molecules has been long neglected. We have 



come across only three previous, very brief suggestions 



experiments now underway at Texas A&M University 111 Oh . 



[|l2l - ll4|1 . Our treatment accompanies 
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In Sec. II, in order to assess the major role of velocity spreads in merged beams, we evalu- 
ate the average relative kinetic energy, (E R ) and its rms spread AE R by integrating over velocity 
distributions familiar for molecular beams. Reduced variable plots are provided that display the 
dependence on the ratio of most probable velocities in the merged beams, their velocity spreads, 
and merging angle. In Sees. Ill and IV we discuss experimental prospects, limitations, and options. 



H. AVERAGES OVER VELOCITY DISTRIBUTIONS 

For beams with lab speeds Vi and V 2 intersecting at an angle 6, the relative kinetic energy is 

E R = x -ix iy\ + v 2 2 -2v { v 2 cose) ( 1 ) 

with fi = m\m%l{m\ + m^) the reduced mass. For merged beams it is feasible to restrict the angle 
6 to a small spread, fixed by geometry and typically only about a degree or so about 9 = 0°. Here 
we evaluate the average of E R over the beam velocity distributions, 

(E R ) = ^ [<V 2 > + <V 2 2 ) - 2<V 1 ><V 2 >cos0] (2) 

and the rms spread, 

AE R = [(E 2 R )-(E R ) 2 ] m . (3) 

These require only (V k ) , with k = 1 - 4, for the individual beams. We obtain analytic expres- 
sions for averages over velocity distributions for beams formed by effusive flow, by supersonic 
expansion, and by a rotating supersonic nozzle. Figure 1 illustrates these distributions. For the 
supersonic beams, three widths are shown; the broadest (~ 10%) is typical, the narrowest (~ 1%) 
is near the best achieved in Stark or Zeeman molecular decelerators exploiting phase stability and 
transverse focusing QVL eD- Results given here pertain to molecular flux distributions; as noted in 
an Appendix, they are readily adapted for number density distributions. 

A. Effusive beams. 

The flux distribution, 

F(V) = V 3 exp[-(V/of) 2 ] (4) 
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is governed by a single parameter, a = (2k B T /m) 1/2 , with k B the Boltzmann constant, and T the 
source temperature. The averaged powers of the velocity are 

Mt)dtl F(t)dt = T[(k + 4)/2] (5) 
o Jo 

with t = V/a and T(z) the Gamma function. The most probable velocity, V mp /a = (3/2) 1/2 = 
1.224, and the rms velocity spread is (AV/a) rms = [2 - (|tt 1/2 ) 2 ] 1/2 = 0.483. The average beam 
kinetic energy is 

(E BK )/(\ma 2 ) = ((V/a) 2 ) = 2 (6) 
and the rms kinetic energy spread is 

AE rms /(\ma 2 ) = [((V/a) 4 ) - ((V/a) 2 ) 2 ] 112 = [Y(4) - T 2 (3)] 1/2 = [6 - 4] 1/2 = 1.41 (7) 

The spread thus is comparable to the beam kinetic energy. 

For a merged pair of effusive beams, it is convenient to define a\ = ancosc/), a2 = a^sin^, 
with a u = (a\ + ff 2 ) 1/2 , and (v k ) = ((V 1 /a 1 ) k ) = ((V 2 /a 2 ) k ) = T[(k + 4)/2]. The relative kinetic 
energy is 

(E R )/(\fxa 2 n ) = (v 2 ) - 2cs(v) 2 cos0 (8) 
with c = cos0 and s = sin0. The rms spread is 

AE R /(^aj 2 ) = [A- Bcos9 + Ccos 2 ^] 1/2 (9) 

with 

A = (c 2 + s 2 ) (<v 4 ) - <v 2 ) 2 ) 
B = 4cs((v)(v 3 )-(v) 2 (v 2 )) 
C = 4cV((v 2 ) 2 -(v) 4 ) 

Figure 2 shows how (E R ) and AE R vary with the ratio of beam velocities, which is proportional 
to a 2 /a\ = tan0. The curves given are for intersection angles near zero, pertinent for merged 
beams. For matched beam velocities, with a 2 = a\, both (E R ) and AE R are smallest. There, 
for small 6, (E R ) and AE R are less than the nominal relative kinetic energy, E° R = j/u(a 2 + a^), 
for a perpendicular collision (6 = 90°) by factors of only about 5 and 3, respectively. As a 2 
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decreases below ct\, (E R ) and AE R increase in nearly parallel fashion. Accordingly, whether or not 
the merged beam velocities are closely matched, AE R exceeds (E R ) appreciably (by ~ 40% when 
a 2 = a\). The resolution of the relative kinetic energy hence is worse than seen in Eqs. (6) and (7) 
for the beam kinetic energy itself. 

The velocity distribution for effusive beams is the same as for a bulk gas. Thus, (E R ) and AE R 
for gas mixtures cooled by cryogenic means can be obtained from the merged beam results by 
merely setting cos# = 0, equivalent to integrating over all angles of collision. 

B. Stationary supersonic beams. 

A standard approximation, 

F(V) = V 3 exp{-[(V - u)/Av] 2 } (10) 

for supersonic beams characterizes the velocity distribution by the flow velocity u along the cen- 
terline of the beam, and a width parameter Av = a\\ = (2k B T\\lm) xl2 , where T\\, termed the parallel 



or lon gitu dinal temperature, pertains to the molecular trans 



ational motion relative to the flow ve- 



locity [|15|1 . According to the thermal conduction model um , T\\IT$ is determined by the pressure 
within the source, P , the nozzle diameter, d, and the heat capacity ratio, y = C P /C V - Likewise, 
the flow velocity is given by 

u = (2k B T Q /m) l/2 [y/(y - 1)] 1/2 [1 - (7y7/ )] 1/2 (11) 

Analytic results for the velocity averages are readily obtained, 

((V/uf) = P k+3 (x)/P 3 (x) (12) 

with x = Av/u the ratio of velocity width to flow velocity, and 

y~»00 

P n (x)= fcxp{-[(t-l)/x] 2 }dt (13) 
Jo 

with t = V/ u. In the Appendix we give exact analytic formulas for the P„(x) functions; Table I 
provides polynomial approximations; for x < 0.3 these are accurate to better than 0.03%. 

The most probable velocity is V mp = ^u{l + [1 + 4x 2 ] l/2 } ~ u{\ + x 2 ). The rms velocity spread 
is (AV/u) rms = [P5/P3 - (Pa/P3) 2 ] 1/2 - The corresponding average beam kinetic energy is 

(E BK )/amu 2 ) = P 5 /P 3 = ((V/u) 2 ) (14) 



and the rms spread in kinetic energy is 

AE BK /(\mu 2 ) = [P 7 /P 3 - (P 5 /P 3 ) 2 ] m (15) 
For a pair of merged supersonic beams, the averaged relative kinetic energy is 

(E R )l(\fiu\ 2 ) = c 2 (v 2 ) + s 2 {v\) - 2cs<v 1 > <v 2 )cos£ (16) 
where u\ 2 = (u\ + u 2 2 ) 112 . The rms spread is 

AE R / (l/uu 2 2 ) = [A - Bcos9 + Ccos 2 #] 1/2 (17) 

with 

A = c 4 [<vl> - (v 2 ) 2 ] + / \{v\) - (v 2 2 ) 2 ] 

B = 4cs [c 2 (<v?> - < Vl ><v?>) <v 2 ) + s 2 ((vl) - (v 2 )(v 2 2 )) < Vl >] 

C = 4cV [(vlXvb-ivtfivrf] 

These expressions are akin to Eqs. (8) and (9), with a\ and a 2 replaced by u\ and u 2 . Here 
(vj ) = ((Vi/Ui) k ), given by Eq.(12), may differ for the two beams (i = 1, 2) if their velocity widths 
(Xj = AVi/ui) differ. 

Figure 3 displays the dependence of (E R ) and AE R on the ratio of flow velocities of the beams, 
W2/W1 = tan0, and their velocity widths. Curves are shown for 8=1° and 2°, to illustrate that 
the dependence on the merging angle is weak if the velocity widths of the beams are fairly large 
(cf. Fig. 2) but for (E R ) becomes significant if the velocity widths are small (cf. Fig. 6 below). 
The dependence on the magnitude of the flow velocities is included simply by adopting units for 
(Er) and aEr that compare them with the relative kinetic energy for collisions at about the most 
probable beam velocities (nearly equal to u\, u 2 ) at right angles {6 = 90°), or equivalently in 
a bulk gas. We designate that by E° R . As seen in panel (a), if the merged beam velocities are 
precisely matched {u 2 ) the collision kinetic energy (E R ) is very sensitive to the 

velocity widths. In contrast, when the beam flow velocities are unmatched by more than about 
15% (i.e., u 2 lu\ < 0.85), the ratio {E R )/E° R becomes nearly independent of the velocity widths 
and grows larger as the unmatch increases. Panel (b) shows that regardless of whether the beam 
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velocities are matched or not, the spread in relative kinetic energy, AE R , varies strongly with the 
velocity widths. Panel (c) plots the ratio AE R /(E R ), which defines the energy resolution. Whereas 
for closely matched beams (E R ) is minimal, AE R then approaches its maximal value. Indeed, for 
matched beams with x > 0.05, the resolution ratio, AE R /(E R ), is near 2 1/2 ; that is just as poor as 
found in Fig. 2 for effusive beams. 

To improve the resolution ratio for matched beams requires narrowing the velocity widths. 
Table II compares, both for the single beam and matched merged beams, effects of reducing the 
spread from x = 0.1 to 0.01. For the single beam, the change in average kinetic energy is very 
slight, whereas the rms spread in the kinetic energy shrinks tenfold. For the merged beams, the 
relative kinetic energy is lowered by a factor of 25, and its rms spread by a factor of 100; so 
the resolution ratio is only improved fourfold. The resolution ratio can be lowered further if the 
merged beam velocities are unmatched, but that raises the averaged relative kinetic energy. Figure 
4 displays the trade-offs involved. To obtain optimally low (E R ) requires nearly exact matching; 
that can provide (E R )/E° R = 2 x 10~ 4 for x = 0.01 or 3 x \0 4 for x = 0.02. However, for exact 
matching, the resolution ratio is only AE R /(E R ) = 0.35 for x = 0.01 and surges to 0.8 for x = 0.02. 
To attain resolution of 0.2 or 0.3 even with x = 0.01 requires u 2 /ui = 0.91 or 0.94 and hence would 
increase (E R )/E° R to ~ 5 x 10 -3 or 2 x 10~ 3 , respectively. The upshot is, to improve the resolution 
ratio by a factor of less than 2 (from 0.35 to 0.2) by unmatching, requires increasing (E R ) by a 
factor of 25. 
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TABLE I: Approximate integrals for velocity averages/' 



n A„(x,y) 

2 4(2y 2 + x 2 )x 

3 4(2/ + 3x 2 )xy 

4 2(4/ + 12x 2 y 2 + 3x 4 )x 

5 2(4/ + 20x 2 y 2 + I5x 4 )xy 

6 (8/ + 60x 2 / + 90x 4 y 2 + \5x 6 )x 

1 (8y 6 + 84x 2 / + 2\0x 4 y 2 + I05x 6 )xy 



"The P„ functions: P„(x) = P„(x, 1); P„{x,y); P n (x,y,z) denned in Eq.(13), and in Eqs.(Al) and (A3) of the 
Appendix, respectively, are all well approximated by P„ = ( -0r/8)A„(jr, y). For x < 0.3 and y > 0.5, the error in 
this approximation for P„ is < 0.03%. 



TABLE II: Comparison of averages and spreads." 



Quantity 


Formula 


x = 0A0 


x - 0.01 


V m p/u 


~ (1 + x 2 ) 


1.010 


1.0001 


AV/u 


[(P 5 /P3)-(P 4 /P3) 2 ] 1/2 


0.070 


0.0071 


AV/V mp 




0.070 


0.0071 


{E BK )l{\mu 2 ) 


P5/P3 


1.035 


1.0003 


AE BK /(\mu 2 ) 




0.143 


0.0141 


AE BK /(E BK ) 




0.14 


0.014 


{Er)K\hu 2 2 ) 


Eq.(16) 


0.0051 


0.00020 


AE R l{\nu 2 2 ) 


Eq.(17) 


0.00697 


0.000071 


AE R /(E R ) 




1.37 


0.35 



"For supersonic beams with un = u\ and 6 = 1°; see Eqs. 14-17 and Figs. 3, 4, and 6. 

C. Rotating supersonic beams. 

For a supersonic beam from a rotating source, the velocity distribution, 

F(V) = V\V - V rot )cxp{-[(V - w)/Av] 2 } (18) 
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involves the peripheral velocity of the source, V rut , which enters both in the V - V rot factor and in 



12D. 



the flow velocity in the laboratory frame, w = u + V rot , with u again relative to the nozzle 113 
In the slowing mode, when the rotor spins contrary to the beam exit flow, V rot < 0; in the speeding 
mode, V > V rot > 0. Also, in integrating over the velocity distribution of Eq. (18), the lower limit 
depends on the sign of V rot . In the slowing mode, the lower limit is small (typically a few m/s) 
and represents a minimum, V > V swat , necessary to allow molecules to escape swatting by the 
rotor. As noted in the Appendix, the swatting correction is extremely small, so we take V swat = 0. 
In the speeding mode, the lower limit may be large, since V > V rot is required. Here, in addition 
to x = Av/u, it is convenient to specify two additional variables: y = w/u and z = V rot /u. The 
Appendix gives analytic results for the integrals, which are denoted by P n (x, y) for the slowing 
mode and P n (x, y, z) for the speeding mode. The corresponding velocity averages are given by 

((V/u) k ) = (P k+3 - zP k+2 )/(P 3 - zP 2 ) (19) 

where P n = P n {x,y) or P n (x, y,z) for slowing or speeding (z < or z > 0), respectively. 

For merged beam experiments with neutral atoms or molecules, pairing a stationary supersonic 
source with a rotating source facilitates adjusting the relative flow velocity. For a stationary source, 
to adjust the flow velocity requires changing the temperature or, if the beam species of interest is 
seeded in a carrier gas, changing the seed ratio. That is awkward and imprecise. For the rotating 
source, the lab flow velocity, w, can be scanned easily and precisely by merely changing the rotor 
speed. The velocity width for a beam from the rotating source is likely to be wider, especially in 
the slowing mode (as Av/w > Av/u). However, it will often be desirable to further narrow the 
velocity spreads in both beams by means of electric or magnetic fields. Such narrowing operations 
can benefit from having the peak intensities of the beams preselcted to occur near the desired 
relative velocity. 

The merged beam (E R ) and AE R can be obtained from Eqs.(16) and (17), using (v^) from 
Eq.(12) and (v^) from Eq.(19). Figure 5 shows, for both slowing (y < 1) and speeding modes 
(y > 1), how results compare with those in Fig. 3 for a pair of stationary supersonic beams (y = 
1). In panels (a) and (b), prominent minima occur where the flow velocities match: w 2 = yu 2 = u\, 
hence u 2 /u\ = tan^> = l/y. In panel (c), the resolution ratio, AE R /(E R ), is least good (largest 



values) at the matching locations, but improves rapidly with modest unmatching, dropping to near 

0. 1 within ~15% of the matching peaks. 

D. Distribution of relative kinetic energy. 

As seen in Figs. 3-5, to attain both low (E R ) and a small resolution ratio AE R /(E R ) requires 
that the merged beams have narrow velocity spreads. Even for x ~ 0.01, to get a resolution ratio 
below ~0.35 requires that the most probable velocities differ somewhat. To examine further the 
competing aspects, we computed numerically the distribution P{E R ) of relative kinetic energy and 
its variation with the velocity spread, extent of unmatching, and the merging angle. 

Figure 6 shows results for merged supersonic beams, both from stationary sources; results with 
one beam from a rotating source are similar. Panel (a) shows P(E R ) distributions for beams having 
velocity spreads of only 0.01 and merging angle 6=1° with ratios of flow velocities u%/u\ = 

1, 0.95, 0.90, and 0.85. Most striking is the form of P(E R ) when the beam velocities are both 
narrow and matched (x = 0.01, «2 = Then the lower limit of the relative kinetic energy, 
E R /E° R < 1.5 x 10^ 4 , is sharply defined and P(E R ) resembles qualitatively a Poisson distribution. 
When the velocities become more and more unmatched P{E R ) becomes approximately Gaussian. 
In contrast, when the beam velocities are broader (e.g., x = 0.10), the P(E R ) distribution is roughly 
Gaussian for both matched and modestly unmatched beam velocities. 

Recasting Eq.(l) in the equivalent form, 



makes evident why the lower limit for E R is sharply defined for merged beams with narrow, 
matched velocities. Then the first term in E R becomes very small, yet as long as 6 + 0, the 
second term will be appreciable because the parent beam velocities are large. That second term 
will also be much less sensitive to velocity spreads; in Fig. 6, panels (b) and (c) exhibit its role. 

III. EXPERIMENTAL PROSPECTS 

In study of chemical reactions at cold (< 1 K) energies under single-collision conditions, con- 
straints enter that differ from "warm" (typically > 300 K) molecular beam experiments. For reac- 
tions without electronic or vibrational excitation, cold reactive collisions cannot occur unless an 




(20) 
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activation barrier is absent or exceptionally low or thin. Such reactions are usually exoergic. The 
deBroglie wavelengths, long for the reactants, thus become very short for the products: "waves 
in, particles out" 11180 . Hence, disposal of energy and angular momentum among product states is 
virtually the same as in the warm regime. Also, as s-waves predominate in the entrance channel, 
the product angular distribution for cold collisions is typically isotropic or nearly so. Accordingly, 
dynamical features observable in warm collisions are not revealed in cold collisions. In the cold 
chemistry regime, the primary information to be sought is the integral reaction cross section and 
its dependence on relative kinetic energy or reactant excitation and/or interactions with external 
fields. 



Merged beams are well suited for measurements of total reaction cross sections p 1.]. The ad- 
vantage at very low collision energies becomes immense as compared with experiments that slow 
both reactants, because the centerline flux in a molecular beam (whether effusive or supersonic) is 
proportional to the most probable velocity in the beam. Although merging imposes a very small 
angular spread, that is compensated because the reactant beams meet in a pencil-like volume, com- 
parable in size to that in typical crossed-beam experiments. Merged beams are also congenial for 
use with pulsed sources, either stationary or rotating, which provide much higher intensity than 
continuous sources. Moreover, the option of overlapping reactants or not in time enhances S/N 
discrimination. Intensity is fostered too by the compactness of merged-beam geometry which also 
somewhat simplifies incorporating auxiliary apparatus. 

In order to facilitate obtaining estimates for any choice of reactant atoms or molecules and 
velocity conditions, we have used unitless ratios in Figs. 2-6 and Table II. Here, in illustrative dis- 
cussion of experimental options, we revert to customary units: energy in mK (or K); mass in amu, 
velocity in meters/sec. The merged beam technique can be used in many variants; we will con- 
sider just two that serve to exemplify basic aspects. We refer to these, somewhat whimsically, as 
"ubiquitous" (ubiq) and "utopian" (utop), to contrast what can be done using ordinary supersonic 
beams with that requiring much more ambitious, sophisticated apparatus. 

Table III displays the chief distinction: for ubiq, x ~ 0.1 is typical, whereas for utop, x ~ 0.01 is 
as yet a challenging goal. The corresponding ranges of averaged collision energy and rms energy 
spread attainable are shown for either matched or 10% unmatched merged beams. As the overall 
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energy scale is governed by E° R = \n{u\ + u\), for either ubiq or utop to obtain the lowest (E R ) 
and aE r requires that the reduced mass and most probable beam velocities be as small as feasible. 
Table III pertains to // = 1 amu, available for reactions of an H atom with any considerably heavier 
reactant partner. The range of flow velocities considered, u ~ 300 to 600 m/s, extends as high as 
likely to be used. Even for light reactants, u ~ 300 m/s can be obtained by seeding in xenon carrier 
gas. That reduces intensity by a factor of typically ~ 100, but is routinely done as a precursor to 
most current methods for slowing molecules. Merged beams, especially if pulsed, can much better 
afford such a drop in intensity because further slowing is not required. 



TABLE III: Illustrative (E R ) and AE R values. 



Mi (m/s) 


U2 (m/s) 


E° R (K) 


(E R ) (mK) 
jc = 0.1 jc = 0.01 


aEr (mK) 
jc = 0.1 x = 0.01 


600 


600 


43 


215 


9 


290 


3 


500 


500 


30 


150 


6 


200 


2 


400 


400 


19 


95 


4 


130 


1.3 


300 


300 


11 


55 


2 


75 


0.7 


600 


540 


39 


390 


195 


460 


39 


500 


450 


27 


270 


135 


320 


27 


400 


360 


17 


170 


85 


200 


17 


300 


270 


10 


100 


50 


120 


10 



fl Kinetic energies pertain to reduced mass fi — 1 amu and beam merging angle 6-1°, with veloc- 
ity spreads (xi = X2) the same in both beams. Entries derived from Table II, but rounded to one or two 
digits. Energy (mK) = 0.0605 mass(amu) [Velocity (m/s)] 2 . Corresponding deBroglie wavelength /l«(nm) = 
395[mass(amu)E(mK)]- 1/2 . 



As seen in Table III, the ubiq mode is capable of reaching (E R ) as low as ~ 55 mK for H 
atom reactions. But, as evident in Figs 3 and 5, the rms energy spread is larger, so the resolution 
is miserable. In the cold collision regime, however, this is not as severe a limitation as might 



be expected. According to Wigner's threshold law II 1911 . for any exoergic inelastic or reactive 
collision at sufficiently low energy the cross section becomes proportional to l/V R , so the rate 
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coefficient becomes independent of the relative kinetic energy. State-of-the art quantum scattering 
calculations [20] for a variety of inelastic and reactive processes indicate the Wigner regime is 
typically attained when E R drops below 100 mK (< 10~ 5 eV), which corresponds to a deBroglie 
wavelength of ~ 40 nm or longer for H atoms. 

To achieve the utop mode by shrinking the velocity spreads to only 1% of the most probable 
velocities, yet preserve adequate intensity, the best means presently in prospect appears to be 
deceleration by multistage Stark or Zeeman fields fl9|], likely exploiting transverse focusing and 



14J] . Table III presumes the resultant 



phase stability together with synchrotron-style storage rings 
velocity distributions resemble Eq.(lO). If so, the utop mode should be able to lower (E R ) to ~2 
mK for H atom reactions, while reducing the energy spread to 0.7 mK, thereby improving the 
resolution to 35%, a rather modest level. As seen in Fig. 4 and Table II, in order to improve to 
the resolution to 20% would require that the beam velocities be unmatched by ~10%, but that 
pushes (Er) up to 50 mK. As the Wigner regime is readily accessible for H atom reactions, even 
for ubig, the elaborate experimental effort to attain utop would not be justified. Obtaining decent 
resolution would become much more worthwhile for reactions with considerably larger reduced 
mass, however. Then, even for utop, the Wigner regime would lie mostly beyond the accessible 
range. 



IV. DISCUSSION 

In summary, the merged-beam technique has several inviting virtues. Foremost is the capabil- 
ity to study cold collisions with warm beams. In contrast to methods that require slowing both 
collision partners, the beam intensities can be much higher and the variety of molecular species 
used much wider. For collisions with small reduced mass, especially H or D atoms with heavy 
molecules, relative kinetic energy below 100 mK can be attained, using mostly standard, fairly 
simple molecular beam apparatus. However, to get low kinetic energy for collisions with large 
reduced mass, and/or to obtain decent energy resolution, will require making the beam velocity 
spreads very small (~1%), a challenging task which entails developing far more elaborate appara- 
tus. 
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The merged-beam experiments implemented in our laboratory aim both to test the capabilities 
of ordinary, basic apparatus JlO| and to explore in the cold regime H atom reactions that have 
been well studied in the warm domain. We note some indicative aspects. Pulsed supersonic 
sources are operated at high input pressures to enhance the familiar features of supersonic beams: 
high intensity, narrowed velocity spreads, drastic cooling of vibration and rotation. One source is 
stationary, the other rotating to enable readily adjusting its beam velocity to match that from its 
sedate partner. The rotor source is suitable for any fairly volatile molecule, while the stationary 
source can be used to generate species that must be produced from precursors, such as hydrogen, 
oxygen, or halogen atoms or free radicals. Initial experiments use the H + N0 2 — > OH + NO 
reaction. The H (or D) beam is supplied by dissociation of H 2 (or D 2 ) in an RF discharge source 



of exceptional efficiency, evolved from a design used in low-temperature NMR [|2 The H beam, 
seeded in Xe or Kr, has flow velocity in the range of 350 - 450 m/s and spread ~10%. The rotating 
source provides the N0 2 beam, currently with estimated velocity spread ~20%. With the beams 
merged at 420 m/s and 6 = 1.5°, the predicted (E R ) ~ 210 mK and AE R ~ 250 mK. That is well 
within the "cold" realm (< 1 K), whereas the kinetic energy is — 11 K for the H beam and ~490 K 
for the N0 2 beam. The corresponding collisional deBroglie wavelength A R ~ 27 nm. 

Here we have considered only effusive and supersonic beams, the most familiar. There is, 
however, a growing repertoire of other types, developed with focus on slowing. Some may offer 
properties advantageous for merged-beams, which shift the focus to narrowing velocity spreads. In 



particular, we mention the recently developed cryogenically cooled buffer gas beams [|22ll . which 
operate in the intermediate regime between effusive and supersonic flow. There the Knudsen num- 
ber (essentially the ratio of Reynolds number and Mach number) is Kn ~ 1 - 10 2 , whereas for 
effusive beams Kn > 1 and for supersonic beams typically Kn < 10" 3 . In a prototype case pro- 



ducing a ThO beam [23], a continuous gas flow maintains a stagnation density of Ne buffer gas in 
thermal equilibrium within a cell cooled to 18 K. A solid Th0 2 target is mounted to the cell wall. A 
YAG laser pulse vaporizes part of the target and ejects ThO molecules that are cooled by collision 
with the Ne buffer as they flow out of the cell. A key effect is "hydrodynamic entrainment" that 
occurs when the time for the ThO molecules to exit the cell is less than the diffusion time to the 
cell walls. Remarkably, the molecular beam has intensity exceeding that for a supersonic beam, 
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very low internal temperature, and its most probable velocity, governed by the cooling Ne, is < 
200 m/s. The velocity spread is broad, but the high intensity and low speed would permit much 
narrowing via velocity selection. Such cryogenically cooled and buffered beams seem likely to 
become widely applicable, perhaps along with merged-beams. 
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APPENDIX: INTEGRALS FOR VELOCITY AVERAGES 



The averages, (V k ), over the velocity flux distributions specified in Eqs. (10) and (18) for 
supersonic beam sources involve two related integrals 

v21 



XOO ( t \ 

fcxp -(-^) 



dt 



r— 2 

[1 + Erf(y/x)] A n (x,y) + ^-exp [-(y/xf] B n (x,y) 



(Al) 
(A2) 



P n (x,y,z) 



f"exp 



t-y 



dt 



(A3) 



I— 2 

= [1 + Erf(l fx)] A n (x, y) + ^-exp [-(1 1 xf] [B n (x, y) + C n (x, y, z)] (A4) 



where, n = k + 3, Erf is the error function, and x = Av/u, y = w/u, z = V rot /u = y - 1. For a 
stationary source, only P n (x, 1) enters, as y = 1 and z = 0; for a rotating source, P n (x,y) pertains 
to the slowing mode, with < y < 1 and z < 0, and P„(x,y, z) to the speeding mode, with y > 1 
and z > 0. Analytical formulas for the A n (x,y), B„(x,y) and C n (x,y, z) functions can be obtained 
by evaluating the P n in successive steps using integration by parts. Table I of the text gives, for n 
= 2 to 7, explicit expressions for A n {x, y); Table IV below gives B,,(x, y) and C n (x, y, z). 
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In the range of interest here, x < 0.3 and y > 0.5, the first term in Eqs. (A2) and (A4) is 
much larger than the others, and the error function is very close to unity. Accordingly, good 
approximations for P n (x,y) and P n (x,y,z) can be obtained by dropping the exponential part and 
replacing the Erf function by unity for compensation; thus, 



P n (x,y) « P n (x,y,z) * ^-A n (x,y) = V^J] 



[n/2] 



s=0 



(25-1)!! Jis+1 „_ 2 , 



(A5) 



For x < 0.3 and y > 0.5, the error is < 0.03%. We note that the swatting correction mentioned 
under Eq.(18) of the text may be obtained from AP n = P n {x,y) - P n (x,y,z*) with z* = V swat /u. 
However, in the range of interest, AP n is extremely small, much less than 0.03% so taking V swat = 
is quite justified. 

The velocity averages are given by 



((V/uf) = (P k+3 - zP k+2 )/(P3 ~ zP 2 ) 



(A6) 



a generic formula that includes Eqs. (5), (12), and (19) of the text. For a stationary supersonic 
source, P n = P n (x) = P n (x, 1); for a rotating source, P n = P n {x,y) or P n (x,y,z) for slowing 
or speeding, respectively. The effusive beam case corresponds to x = a, y = 0, z = 0, with 
P n (x,0) = ja n+l Y[(n + l)/2]. Results for averages over number density distributions rather than 
flux distributions can be obtained merely by setting n = k + 2. 

TABLE IV: Formulas for B n (x,y) and C n (x,y, z). 



B n (x,y) 



C n (x,y,z) 



Ay 

4(v 2 + x 2 ) 
2(2y 2 + 5x 2 )y 
2(2y 4 + 9x 2 y 2 + 4x 4 ) 

(4v 4 + 2Sx 2 y 2 + 33x 4 )y 



7 4/ + 40x 2 y 4 + 87/v 2 + 24x 6 



2 
3 
4 
5 



4z 
4z(z + v) 
2[2z 3 + 2vz 2 + (2v 2 + 3x 2 )z] 
2[2z 4 + 2vz 3 + 2(v 2 + x 2 )z 2 + (2y 3 + lx 2 y)z] 
4 Z 5 + 4yz 4 + 2 (5 X 2 + 2;y 2) z 3 + 2 (2 V 3 + 9x 2 v)z 2 
+(4/ + 24x 2 y 2 + 15x 4 )z 

4 Z 6 + 4yz 5 + 4 q x i + y2 )z 4 + 2( 2v 3 + l l x 2 v)z 3 
+2(2v 4 + 15x 2 v 2 + 12x 4 )z 2 + (4y 5 + 36x 2 y 3 + 57x 4 y)z 
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FIG. 1 : (Color online) Velocity distributions of molecular flux, F(V), for beams formed (a) by effusive flow, 
(b) supersonic expansion from a stationary source, and (c) a rotating supersonic source, defined by Eqs.(4), 
(10), and (18), respectively. Parameters for (a) are Tq = 300^; a = 250 m/s; for (b) and (c) flow velocities 
are u = 400 m/s and w - u + V rot = 250 m/s and widths Av/u = 0.01, 0.05, 0.10. For (b) and (c) the widths 
also influence somewhat the most probable velocity. 
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FIG. 2: (Color online) Properties of merged effusive beams, from Eqs.(8) and (9), obtained by averaging 
over beam velocity distributions, (a) Averaged relative kinetic energy; (b) its rms spread; (c) ratio of spread 
to averaged kinetic energy, defining the available resolution. Ordinate energy scale is E° R - \[i{a\ + or|); 
abscissa scale pertains to ratio a^Mi ranging from 0-2 = \a\ to = a\. Results are shown for intersection 
angles of 9 - 1°, 2° and 5°. 
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FIG. 3: (Color online) Properties of merged supersonic beams, from Eqs.(16) and (17), obtained by aver- 
aging over beam velocity distributions, (a) Averaged relative kinetic energy; (b) its rms spread; (c) ratio 
of spread to averaged kinetic energy. For (a) and (b) the ordinate energy scale is E° R = \[i(u\ + w|). The 
abscissa scale pertains to the ratio uilu\, which ranges from ui = to ui - u\. Results are shown for 
an intersection angle of 6 = 1° and four sets of velocity spreads: x\ = X2 = 0.01; 0.05; 0.1; and x\ =0.1, 
%2 = 0.2. Dashed curves included for the x = 0.01 case are for 8 = 2°. 
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FIG. 4: (Color online) Variation with extent of matching of merged beams flow velocities, uilu\, for relative 
kinetic energy, (Er)/E° r (left ordinate scale) and resolution ratio, aEr/{Er) (right ordinate scale). Curves 
pertain to velocity spreads in beams of x\ - x-i - 0.01 and 0.02 and merging angle of = 1°. 
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FIG. 5: (Color online) Properties of merged supersonic beams, one from a stationary source, the other from 
a rotating source. As in Fig. 3, (a) shows the averaged relative kinetic energy; (b) its rms spread; (c) the 
resolution ratio, spread to averaged kinetic energy. Results shown are for 9=1° and x - 0.01. Both (a) and 
(b) exhibit pronounced minima where the matching condition holds: W2 = u\ and hence U2/u-[ = l/y. Full 
curves show results for slowing mode, with y — 4/5, 2/3, and 1/2; dashed curves are for speeding mode, 
with y = 5/4, 3/2, and 2. 
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FIG. 6: (Color online) Distributions of the relative kinetic energy, P(Er), for merged supersonic beams 
from stationary sources. Values of (Eg) are indicated by black dots. Abscissa scale is in units of E° R = 
\n{u\ + ^2)' as use d i n Figs. 3-5. (a) For merging angle 6-1° and various ratios of the flow velocities, 
ui/ui = 1 to 0.85 and velocity widths xi = x% = 0.01. Note log-log plot is used, (b) For matched flow 
velocities, U2/U1 = 1 and x = 0.01, but merging angle varied to illustrate its role in Eq.(20). (c) For flow 
velocities unmatched by ~15% (u\U2 = 0.85) and x = 0.01 with 8 = 1°, 2°, 3°, to illustrate the reduced role 
of the merging angle when the first term in Eq.(20) becomes predominant. 
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